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Symmetry restoration in a theory of a self-interacting charged scalar field at finite temperature 
and in the presence of an external magnetic field is examined. The effective potential is evaluated 
nonperturbatively in the context of the optimized perturbation theory method. It is explicitly 
shown that in all ranges of the magnetic field, from weak to large fields, the phase transition is 
second order and that the critical temperature increases with the magnetic field. In addition, we 
present an efficient way to deal with the sum over the Landau levels, which is of interest especially 
in the case of working with weak magnetic fields. 
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Phase transition phenomena in spontaneously broken quantum field theories have long been a subject of importance 
' (-H ' , and interest due to their wide range of possible applications, going from low energy phenomena in condensed matter 
O i' systems to high energy phase transitions in particle physics and cosmology (for reviews, see for example ^hSI])- 
' \ In addition to thermal effects, phase transition phenomena are also known to be triggered by other external effects, 
like, for example, by external fields. In particular, those changes caused by external magnetic fields have attracted 
considerable attention in the past ^] and received reinvigorated interest recently, mostly because of the physics 
I— I. associated with heavy-ion collision experiments. In heavy-ion collisions, it is supposed that large magnetic fields can 
be generated, and the study of their effects in the hadronic phase transition then became subject of intense interest (see 
^! . e.g. for a recent review). Magnetic fields can lead in particular to important changes in the chiral/deconfinement 
QQ ' transition in quantum chromodynamics (QCD) [6] and even the possibility of generating new phases Q 

\ As far the influence of external magnetic fields and thermal effects on phase transformations are concerned, one well 
. known example that comes to our mind is the physics associated with superconductivity, in particular in the context 
■^^^ ' of the Ginzburg-Landau theory ^ . Let us recall in that case thermal effects alone tend to produce a phase transition 
QQ . at a critical temperature where superconductivity is destroyed and the system goes to a normal ordered state. The 
' phase transition in this case is second order. However, in the presence of an external magnetic field, but below some 
T-H , critical value, by increasing the temperature the system undergoes a first order phase transition instead. This simple 
• example already shows that magnetic fields may have influence on the phase transition other than we would expect 
^ from thermal effects alone. There are also other examples of more complex systems where external magnetic flelds 
• '-j may have a drastic effect on the symmetry behavior. Among these effects, besides the possibility of changing the 
order of the phase transition, as in the Ginzburg-Landau superconductor, it can in some circumstances strengthen 
^ ,; the order of the phase transition, like in the electroweak phase transition in the presence of external flelds Q , or there 
- - ' can also be dynamical effects, like delaying the phase transition [1^]. External magnetic flelds alone can also lead to 
dynamical symmetry breaking (magnetic catalysis) [ll| (for an earlier account, see Ref. [13 )■ 

Likewise, nonperturbative effects may affect the symmetry properties of a system, once the external parameters are 
changed, in a way different than seeing through a purely perturbative calculation, or by a mean-field leading order 
description. This is because perturbation theory is typically beset by problems, for example around critical points, due 
to infrared divergences, or at high temperatures, when powers of coupling constants can become surmounted by powers 
of the temperature (see e.g. the textbooks lH, Hi! for extensive discussions). Thus, high temperature field theories 
and the study of phase transitions in general require the use of nonperturbative methods, through which large classes 
of terms need to be resummed . Familiar techniques used to perform these resummations include, for example, ring 
diagram (or daisy and superdaisy) schemes [TBI . [l6j , composite operator methods jl7| and field propagator dressing 
methods [isl . [l9| . Other methods used include also numerical lattice studies and expansions in parameters not related 
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to a coupling constant, like the 1/N expansion and the e-expansion [20|, the use of two-particle irreducible (2PI) 
effective actions [2l|, hard-thermal-loop resummation j23|, variational methods, like the screened perturbation 
theory [l^, [l^ and the optimized perturbation theory (OPT) [2^. Of course, any resummation technique must be 
implemented with care so to avoid possible overcounting of terms and lack of self-consistency. Failure in not following 
this basic care can lead to a variety of problems, like predicting nonexistent phenomena or producing a different 
order for the phase transition. One classical example of this was the earlier implementations of daisy and superdaisy 
schemes, that at some point were giving wrong results, e.g. predicting a first order transition [27] for the theory, 
an unexpected result since the model belongs to the universality class of the Ising model, which is second order. 
These methods have also initially predicted a stronger first order phase transition in the electroweak standard model, 
a result soon proved to be misleading [2^ . These wrong results were all because of the wrong implementation of the 
ring-diagram summation at the level of the effective potential, as clearly explained in the first reference in j28j 

In this work we will analyze the phase transition for a self-interacting complex scalar field model and determine how 
an external magnetic field, combined with thermal effects, affects the transition. All calculations will be performed in 
the context of the OPT nonperturbative method. Our reasons for revisiting here the phase transition in this model 
are two-fold. First because this same model has been studied recently in the context of the ring-diagram resummation 
method (291] . where it was found that the ring-diagrams render the phase transition first order and that the effect 
of magnetic fields was to strengthen the order of the transition and also to lower the critical temperature for the 
onset of the (first order) phase transition. So in this work we want to reevaluate these findings in the context of 
the OPT method. We recall, from the discussion of the previous paragraph, that the ring diagram method requires 
special attention in its implementation and that previous works have already concluded erroneously about its effects 
on the transition. The OPT method has a long history of successful applications (for a far from complete list of 
previous works and applications see e.g. Refs. [s^, [HI and references therein). The OPT method automatically 
resums large classes of terms in a self-consistent way so to avoid possible dangerous overcounting of diagrams. In 
our implementation of the OPT here, we will see that already at the first order in the OPT it is equivalent to the 
daisy and superdaisy schemes. The OPT then provides a safe comparison with these other nonperturbative schemes. 
In particular, to our knowledge, this is the first study of the OPT method when considering the inclusion of an 
external magnetic field. Finally, we also want to properly treat the effects of small magnetic fields (which requires 
summing over very large Landau-levels) in an efficient way, particularly suitable for numerical work. This way we can 
evaluate in a precise way the effects of external magnetic fields ranging from very small to very large field intensities 
(in which case, in general, just a few Landau-levels suffice to be considered). For this study we will make use of the 
Euler-Maclaurin formula and fully investigate the validity of its use as an approximation for the Landau level sums 
for different ranges for the magnetic field. 

The remaining of this paper is organized as follows. In Sec. II we introduce the model and explain the application 
of the OPT method for the problem we study in this paper. It is shown explicitly which terms are resummed by the 
OPT. We also verify that the Goldstone theorem is fulfilled in the OPT. In Sec. Ill we study the phase transition 
in the spontaneously broken self-interacting quartic complex scalar field in the OPT method, first by incorporating 
only thermal effects and then by including both temperature and an external magnetic field. In Sec. IV we study the 
advantages of using the Euler-Maclaurin formula for the sum over the Landau levels. The accuracy and convergence 
of the method is fully examined. We determine that the sum over the landau levels, in the regime of low magnetic 
fields, where typically we must sum over very large levels so to reach good accuracy, can be efficiently performed 
within the very few first terms in the Euler-Maclaurin formula. As an application, an analytical formula for the 
critical temperature for phase restoration, as a function of the magnetic field, Tc{B), is derived. Our final conclusions 
are given in Sec. V. Finally, an appendix is included to show some of the details of the renormalization of the model 
in the context of the OPT. 



II. THE COMPLEX SCALAR FIELD MODEL AND THE OPT IMPLEMENTATION 

In our study we will make use of a self-interacting quartic complex scalar field model with a global U(l) symmetry 
and spontaneously symmetry breaking at tree-level in the potential, whose Lagrangian density is of the standard form, 



L = \d,(t>?+m^W-^W , (2.1) 

where m? > for spontaneously symmetry breaking. It is convenient to write the complex scalar field in terms of 
real and imaginary components, (j) — (0i + i(j)2)/V2. In terms of a (real) vacuum expectation value (VEV) for the 
field, {(j)) = (p/V2, we can, without loss of generality, shift the field around its VEV in the (pi direction. 
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(2.2) 



The Feynman propagators for 0i and (f>2, in terms of the VEV then reads, 



D^, (P) = ^ , , (2.3) 

Using the tree-level VEV for the field, ipQ — 6m^/A in Eqs. p.3p and ()2.4p . 0i is then associated with the massive 
Higgs mode (with mass squared = 2m^ at the tree level) and (f>2 is the Goldstone mode of the field, remaining 
massless throughout the symmetry broken phase. 



A. Implementing the Optimized Perturbation Theory 



The implementation of the OPT in the Lagrangian density is the standard one [30, l31[ , where it is implemented 
through an interpolation procedure. 



^ ^ = E {I (9,^^)" - \^"^" + U'^" - ^ i^f] + , (2.5) 

where f2^ = —m^ + ij^ and ALct.5 is the Lagrangian density part with the renormalization counterterms needed to 
render the theory finite. In Lg, the dimensionless parameter (5 is a bookkeeping parameter used only to keep track 
of the order that the OPT is implemented (it is set to one at the end) and ry is a (mass) parameter determined 
variationally at any given finite order of the OPT. A popular variational criterion used to determine rj is known as 
the principle of minimal sensitivity (PMS), defined by the variational relation (32| 



drj 



= , (2.6) 



which is applied to some physical quantity (J)''^', calculated up to some order-/c in the OPT. The optimum value fj, 
which satisfies Eq. (j2.6p , is a function of the original parameters of the theory and it is in general a nontrivial function 
of the couplings. It is because of the variational principle used that nonperturbative results are generated. Other 
variational criteria can likewise be defined differently, but they produce the same final result for the quantity ^C'), as 
shown recently [33| . 

In terms of the interpolated Lagrangian density, Eq. (|2.5p . the Feynman rules in the OPT method are as follows. 
The interaction vertex is changed from —i\ to —iSX. The quadratic terms 6r]^(l)f/2 in Eq. (|2.5p define new vertex 
insertion terms. The propagators for (jji (i = 1,2) now become 



B. The effective potential in the OPT method 



As in the many other previous applications [31|,|33j, we apply the PMS, Eq. (j2.6p . directly on the effective potential, 
which is the most convenient quantity to study the phase structure of the model. To first order in the OPT, using 
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FIG. 1: Feynman diagrams contributing to the efTective potential to first order in the OPT. Sohd lines stand for the 0i 
propagator, while dashed lines stand for the 02 propagator. A black dot is an insertion of drf . A crossed dot is a mass 
renormalization insertion. 

the previous Feynman rules in the OPT method given above, the effective potential Vcs{ip) is given by the vacuum 
diagrams shown in Fig. [1] 

From the diagrams shown in Fig. [T]and by further expanding the propagators Eqs. (|2.7I) and (|2.8|) in (5, the explicit 
expression for the renormalized effective potential at first order in the OPT (some of the details of the renormalization 
in the OPT method are given in the Appendix A) becomes 



ln(P 
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(2.9) 



where the momentum integrals are expressed in Euclidean space and the regularization is performed in MS scheme, 
where, in the finite temperature only case. 




(2.10) 



where 7e is the Euler-Mascheroni constant, M is an arbitrary mass regularization scale, d = 3 — 2e is the dimension 
of space and a;„ — 2TmT, (n = 0, ±1, • • • ) are the Matsubara frequencies for bosons at temperature T. 

The inclusion of an external magnetic field presents no additional difficulty. For example, without loss of generality, 
we can consider a constant magnetic field in the z-space direction and use a gauge where the external electromagnetic 
field is Afj^ — (0, 0,i?x, 0). The Feynman propagator for a charged scalar field with charge e becomes P^ — — )• 
—uj'^ — E^{pz, B), where Ek{pz, B) is the energy dispersion (for charged scalar bosons) in an external constant magnetic 
field ^], 



Ei^pi + m^ + {2k+l)eB, (2.11) 

where the last term denotes the Landau levels (A; = 0, 1, 2, • • • ). Likewise, the momentum integrals, taking into account 
the degeneracy multiplicity of the Landau levels [35| , are now represented by 

ePi2?^ ^ fe^^M^y f d^-^p. 



i—\'T >^ ^ / ; ■ (2.12) 

k=0 Po=iUJr, ^ / J \ ) 



C. Optimization results and Goldstone theorem 



Let us now apply the PMS condition (12.61) on the effective potential Eq. (|2.9p . We obtain straightforwardly that 
the optimum f\ satisfies the nontrivial (renormalized) gap equation 
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(2.13) 



The extrema of the effective potential are defined as usual, by requiring that the first derivative of the effective 
potential with respect to ip vanishes. This gives us the trivial solution (p = and the solution for the minimum, 
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(2.14) 



We can now verify the effective mass for the field, in particular for the 02 component of the complex scalar field. 
The effective mass for 02 in the OPT is given by 



'■cff,2 
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(2.15) 

where ^2(77) is the (renormalized) field's self-energy, which at first order in the OPT is trivially found to be given by 
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(2.16) 



Using now Eqs. (|2.16p . (|2.13p and (|2.14p in Eq. (|2.15p . we obtain immediately that m^j^ 2=0) which is nothing 
but the result expected due to Goldstone theorem for the symmetry broken complex scalar field model. It can also be 
verified that this result carries out at higher orders in the OPT method, in which case, at some order-fc in the OPT, 
the self-energy is the one at the respective order, S^'^^ entering in the above equations. Previous demonstrations 
of the Goldstone theorem in the OPT were for the linear sigma model [36] and for the SU{2) Nanibu-Jona-Lasinio 
model m. 

Finally, it is noticed from the PMS equation given above, Eq. (|2.13l) . that the OPT naturally resums the leading 
order loop terms of the field's self-energy. So it is quite analogous, at already in the first order in the OPT approx- 
imation, to the ring-diagram resummation [isl [28j . And that this resummation is also performed automatically in 
a self-consistent way, it is clear from the interpolation procedure. While fj enters in all propagators, thus carrying 
the self-energy corrections, the insertions of fj (the term 5rf'4)1/2 in Eq. (|2.5p is treated as an additional interaction 
vertex) subtract spurious contributions at each order, thus avoiding dangerous overcounting of diagrams (this is also 
similar to the procedure used in the screened perturbation theory [2J,[25). Going to higher orders in the OPT resums 
higher order classes of diagrams. For our purposes of comparing with the ring-diagram results, it suffices then to keep 
up to first order in the OPT. 



III. PHASE STRUCTURE AND SYMMETRY RESTORATION AT FINITE TEMPERATURE AND IN 

AN EXTERNAL MAGNETIC FIELD 

We are now ready to study the phase structure and the symmetry restoration in the complex scalar field model at 
finite temperature and in an external magnetic field. It is convenient to first investigate the case of finite temperature 
only, so to later compare with the case when an external magnetic field is added. 

A. Symmetry restoration at finite temperature 

It is convenient to write the explicit expressions for each term entering in the effective potential in Eq. (|2.9p . At 
finite temperature and in the absence of an external magnetic field, using (|2.10p . we have that the first momentum 
integral in (|2.9p becomes (recalling that = —w? + rf') 

-^lHP^-^') = -^l+nT,v), (3.1) 
where we have identified explicitly the divergent term and the finite term, Y{T^ rj), is given by 
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Likewise, for the remaining momentum integrals in (|2.9p we obtain 



where X{T, rj) is given by 
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Next, there is the term coming from the mass countertcrm in Eq. (|2.9I) (see Appendix A), 
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where we have also defined 
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Finally, we have the two-loop contributions shown in Fig. [TJ They all sum up in the OPT expansion to give 



4 



P2 _ Q2 



,x 1 

= — 



3 {4Trf £2 

2A 



.A 



„ vT4/(77)+,5-A2(T,r;) 
3 (47r) 3 



(3.8) 



From Eqs. p.ip . p.3p . p.4p . p.6p and p.Sp . we find the complete expression for the renormalized effective potential 
at first order in the OPT: 



+ Y{T,f^) + 6l^-rf + ^[^^+X{T,f^)]'^X{T,fj) . 



(3.9) 



The phase structure at finite temperature is completely determined by Eqs. p.9p . (|2.13l) and (|2.14p . The optimum 
fj, determined by the PMS criterion, and the vacuum expectation value (VEV) for the scalar field, can now be both, 
respectively, be expressed as 
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f = -v' + ^XiT,fj), 



(3.10) 
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and 



(3.11) 



In the OPT we can exactly compute the critical temperature of phase transition. Using that at the critical point 
the VEV of the field vanishes, <f{Tc) = 0, then from Eq. (I3.10|) we obtain that 



f{T,)^-X{T,,fj{T,)) 



which upon using it in Eq. p. Ill) , we obtain 



= 6—- 4X(r„ fjiTc)) ^m^~ f{T,) = -^''{T,) = 



From the definition of X{T,ri), Eq. p.Sp . we then obtain that 



X{T,,fj{T,)) = ^ dz 



z 
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exp (z) — 1 12 



(3.12) 



(3.13) 



(3.14) 



and when using the above result back in Eq. p.l3p . we obtain the exact result for the critical temperature at this 
first order in the OPT: 



1877?/ 

A 



(3.15) 



This result for is the same predicted before for a scalar field model with two real components [3^, obtained in 
the high temperature one-loop approximation. The result (j3.15p is exact in our OPT approximation and obtained 
independent of any high temperature approximation, as usually assumed in any previous calculations. 
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FIG. 2: Temperature dependence of the minimum of the effective potential ip{T) in the OPT with the PMS optimization 
criterion, ipo is the minimum of the tree-level potential. The parameters used are m/M = 20 and A = 0.00375. 



From Eqs. p.lO|) and (|3.1ip . we obtain directly the behavior of the VEV as a function of the temperature. For 
comparison purposes, we use analogous parameters as adopted in Ref. [2^ ^, where m/M — 20 and A = 0.00375. 



^ Note that in f29il the authors used a different numerical factor for the quartic term in the tree-level potential, which gives rise to the 
extra factor 3/2 in the numerical value for A in our notation. In the approximations used in [29| | , the terms involving the regularization 
scale were not included. Here we express all quantities in terms of the regularization scale M . 
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For these parameters, the authors in Ref. [231 ^^'^ ^ ^^^^ order phase transition. In Fig. [2] we show the resuh for the 
minimum of the effective potential, = ip{T) as a function of the temperature. The effective potential, for different 
temperatures and same parameters as in Fig. [21 is shown in Fig. [3l Note that the effective potential has an imaginary 
part for temperatures below Tc, Eq. p.lSp . since there can be values of (p for which —m^ + becomes negative. 
This happens for values of in between the inflection points of the potential, which defines the spinodal region of 
instability in between the (degenerate) VEV of the field, determined by Eq. p. lip . In Fig. [Sj for the case of T < Tc, 
we show the real part of the effective potential, so to be able to show all the potential, including the spinodal region. 
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FIG. 3: The effective potential (subtracting the vacuum energy at <^ = 0), for the same parameters of Fig. [2] and for 
temperatures: Ti/M = 1380, T2/M = 1390 and at the critical temperature, tJ/M ~ 1385.64. 

It should be noted that the minimum of the effective potential, determined by the coupled Eqs. (|3.10p and (|3.11l) . 
is well defined and unique, where the only solutions for Eq. (jS.llI) are the two degenerate minima for T < Tc, while 
for T > Tc, the only solution is = 0. Therefore, the phase transition can readily be inferred to be second order for 
any set of parameters. It is also clear from the results shown by Figs. [21 and [H that the transition due to thermal 
effects only is of second order. The effective potential does not develop local minima and the VEV of the field varies 
continuously from T = till T = Tc, where it vanishes. This is the expected behavior, since the complex scalar field 
model belongs to the same universality class of a 0{2) Heisenberg model, for which the phase transition is second 
order 

Thus, we see that the OPT method correctly describes the phase transition in the model, capturing the correct 
physics. Furthermore, as an added bonus, we can exactly compute the critical temperature of transition 

B. The phase structure in a constant external magnetic field 

Let us now investigate the phase structure when an external magnetic field is also applied to the system. In the 
presence of a constant magnetic field, the Feynman rules change as shown in Subsec. IIIBI The momentum integrals 
in Eq. (12. 9p are given by (|2.12l) and the dispersion relation for charged particles is given by Eq. (|2.1ip . which takes 
into account the Landau energy levels of a charged particle in a magnetic field^. 

At finite B, the first momentum integral term in (|2.9p becomes: 

HP'-^')-t^[-^) rhA^ ^ ln[.l+pl + n- + [2k + l)eB] . (3.16) 



Note that the dispersion relation Eq. I|2.1ip only applies for the fields in the complex base (0, (/>*), since {<f>i,if>2) are not appropriate 
eigenstates of charge. However, when re-expressing the effective potential in the complex scalar field base, the final result turns out to 
be the same as in Eq. 1 12.911 or 113.91 1. but with the functions X and Y depending now also on the magnetic field. 
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By performing the sum over the Matsubara frequencies in p.l6p . we obtain two terms. One is independent of the 
temperature and the other is for T ^ 0. The T — term is 
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where to write the second hue in p.l7p . we used the analytic continuation of the Hurwitz zeta function (39j . 
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and C'(s, a) is its s-derivative. In writing the final terms in Eq. p.l7p we have also dropped constant terms since they 
do not affect the phase structure. 
The T 7^ part of (PH]) is 
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We can now write Eq. (|3.16p in a form analogous to Eq. p.ip . 
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where Y(T, B,r]) is given by 
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Analogous manipulations leading to Eq. p. 201) allow us to write the remaining momentum integral terms in Eq. 
Uke 
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Actually, Eq. p. 221) is just the derivative of Eq. p.20p with respect to Q'^, as can be easily checked. 
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The final expression for tlie renormalized effective potential at finite temperature and in a constant external magnetic 
field has an analogous form as Eq. p.9p . where by using Eqs. p.20[) and p. 221) in (I2.9p . we obtain 



+ y (T, i?, 77) + ,5 1 -if + ^[p^+ XiT, B, 77)] I XiT, B, r,) . (3.24) 

Finally, the PMS criterion and the minimum of the effective potential are again also of the form as Eqs. p.lOp 
and l|OT|) : 

f = l^'+'^-^X{T,B,n). (3.25) 

and 



2 

m 



^6—-'iXiT,B,fj) . (3.26) 

With Eqs. p.24p . p. 251) and p.26p . we are now in position to study the effects of an external magnetic field in 
the phase structure of the model, in addition to the thermal effects. In Fig. |4] we show the phase diagram for the 
symmetry broken complex scalar field model in the (i?, T) plane. It shows that the magnetic field strengths the 
symmetry broken phase, increasing the critical temperature for phase transition. We must note that this result for 
the phase diagram is very different from that seen in the superconductor or in the scalar quantum electrodynamics 
(QED) case [4^. In contrast to our case, where we are only studying the effects of B and T on the charged scalars, 
in the scalar QED there is also the interactions with the gauge field. In the scalar QED case with a local U{1) gauge 
symmetry, in the symmetry broken phase, because of the Higgs mechanism the gauge field becomes massive and the 
external magnetic field gets screened (the Meissner effect). As the magnetic field is increased, eventually above a 
critical field the phase is restored through a first order phase transition. In the study made in this paper, on the other 
hand, we consider only the case of a broken global U{1) symmetry'^. 

The effect of the magnetic field on the global U (1) symmetry can also be seen in the next two figures. In Fig. [5] we 
have plotted the VEV of the field at fixed T > Tc, thus starting from a symmetry restored phase, 1^ = 0, as a function 
of B. The critical temperature, for the parameters used, is Tc/M ~ 13.42, which fully agrees with Eq. (I3.15|) . and the 
symmetry returns to be broken, (p Q, for a magnetic field above a critical value given by eB^/M'^ ~ 53.93. In Fig. 
inithe effective potential is plotted for the same fixed temperature of Fig. [S]and for values of B below, at and above 
the critical value Be, above which the symmetry becomes broken again (as before, to show the spinodal region of the 
potential, we have plotted only the real part of it). It should be noticed from these results that the phase transition 
is once again second order. The effect of the magnetic field is to enlarge the symmetry broken region by making the 
critical temperature larger, but it does not change the order of the transition. 

Finally, we have again used for comparison a case with analogous parameters as considered in Ref. [29] , and contrary 
to the results found in that reference, we have found again here only a second order phase transition. For the analogous 
parameters used in [2^ (see observations made in previous subsection), m/M = 20, A = 0.00375 and eB/M'^ = 30, 
the change of the critical temperature in relation to the B = Q case shown in Fig. |3]is only to produce a slight shift 
of the critical temperature to a value 0.08% higher. In Fig. [71 we have plotted the VEV of the field for the case 
of zero magnetic field and for a much higher value of the magnetic field, eB/AP — 3000, so to be able to visualize 
the difference. Even so, the difference even for such value of magnetic field is only marginal, changing the critical 
temperature obtained at S = 0, Tc{B = 0)/Af ~ 1385.64, to Tc{B)/M ~ 1396.87 for the value of B used. In ah 
cases we have tested for a nonvanishing magnetic field and other different parameters of the potential, the transition 
is again verified to be second order, with the VEV varying always continuously from its value at T = to zero at 
T = T,(B). 



^ One of the authors (ROR) would like to thank K. G. Klimenko for discussions on this topic and for also pointing him out the possible 
similarity with the symmetry behavior found in this paper due to the magnetic field, with pion condensation seen in the Nambu— Jona- 
Lasinio model with 2 flavor quarks plus baryon and isospin chemical potentials [41II . 
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eB/NT 

FIG. 5: Magnetic field dependence of the VEV of the field, (p{T, B), for a fixed value of T above Tc- The parameters considered 
here are: A = 0.1, m/M = 1 and T/M = 15. 



IV. USING THE EULER-MACLAURIN FORMULA FOR THE SUM OVER LANDAU LEVELS 

When working with quantum field theories in a magnetic field, we need to deal with the sum over the Landau levels. 
At T = this does not present a problem in general, since we can express the expressions in terms of zeta functions. 
We have seen this in the previous section, where all terms at T = were presented in analytical form. However, at 
finite temperature this is not possible in general, therefore, we either need to numerically perform the sums over the 
Landau levels, or find suitable approximations for the expressions. This job becomes easier in the high magnetic field 
regime, with V eB ^ T, for which most of the time suffices to use only the very first terms of the sum, or even just the 
leading Landau level term. Higher order terms are quickly Boltzmann suppressed in this case. But this is not so in 
the low magnetic field regime, with ^/ eB < T , and fi/T < 1, where most of the time requires to sum up to very large 
Landau levels. This has the potential to make any numerical computation to become quickly expensively, specially 
when there are in addition numerical integrations involved, like in our case here (note that we have not made use of 
high temperature approximations, but worked with the complete expressions in terms of the temperature dependent 
integrals). Simple approximate analytical results, which are always desirable to obtain, can also become difficult to 
obtain without a suitable approximation that can be used. One such approximation was discussed in Ref. [1^ , but it 
is valid only for very small magnetic fields. Here we discuss a more natural alternative for dealing with the cases of 
low magnetic fields, based on the Euler-Maclaurin formula. 
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FIG. 6: The effective potential for tlie same parameters of Fig. [5] and for tliree different values of the magnetic field: eB /M^ 
45, eBc/M^ = 53.93 and eB/M^ = 60. 
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FIG. 7: Temperature dependence of the normalized VEV ip{T, B)/^po. The critical temperature for B = is Tf^'^^jM ~ 1385.64. 
Including the effects of the mag netic field, with eB/Af^ = 3000, Tc increase to Tc/M ~ 1396.87. The other parameters used 
here are: A = 0.00375 and m/M = 20. 

The Euler-Maclaurin (EM) formula provides a connection between integrals and sums. It can be used to evaluate 
finite sums and infinite series using integrals, or vice- versa. In its most general form it can be vi^ritten as (4^ : 



1 " /) r 

f{x) dx + i [/(a) + fib)] + E (lyi [f^''^'Hb) ~ f^"'-'Ha) 

B2n+l{{x}) „(2„+l) , . , 

, (2n+l)! ^ ^ '"^ ' 



where bi are the Bernoulli numbers, defined by the generating function 



(4.1) 



exp( 



n=0 



(4.2) 



and Bn{x) are the Bernoulli polynomials, with generating function 
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The notation {x} in B2k+i{{x}) in Eq. (I4.ip means the fractional part of x and f'^^\x) means the fc-th derivative of 
the function. The last term in Eq. (14. ip is known as the remainder. 

Next, we study the reliability of the use of the EM formula as a suitable approximation for the sums over Landau 
levels and we estimate the errors involved in doing so. In the following, we will call the Oth-order in the EM formula 
when only the first term in the RHS of Eq. (|4.ip is kept, the Ist-order when also the second term is kept and so on 
so forth. Thus, for example, the Landau sum part of Y{T, B, 77) in Eq. (|3.2ip . up to second order in the EM formula, 
becomes 



LY 



dz 



■In 



1 - 



^ 7 TO 27r 



E(z,n,T,B,k) 



/•+°° dz ( 










1 - e 


~ J-00 2^1 







-E(z,Q.,T,B,k) 



dk 



eB 



12r2 £:(z,17,T,B,0) [eE{-.^,T,B,o) _ 1] 



(4.4) 



where 



£;(z,f],T,i?,fc) = ^z2 + ^ + (2fc + i)£g . (4.5) 
Likewise, the Landau sum part of X{T, B, rj) in Eq. p.23p . up to second order in the EM formula, becomes 

- _ ^ /■+°° dz 1 1 /■+°° dz f f°° 1 1 

" ^0 i-oo 27r E{z, n, T, B, k) e^^(-.s^^i^^s,/c) _ 1 " 2n [Jo E{z, il, T, B, k) e^^^^^'^'^^k) _ 1 

1 1 



2 E{z,n,T,B,0) [eEi^,^,T,B,o) _ 1] 

eB 1 
12T2 E'^(z,n,T,B,0) [eEi^,n,T,B,Q) _ ^ 



g_E(z,O,T,S,0) 



E{z,n,T,B,0) e^(^AT,i3,o) „ 1 



(4.6) 



In Tab. Uwe assess the reliability of using the EM at each order (we have studied it up to the forth-ordcr) for the 
sums in Eqs. (|4.4p and ()4.6p . for different values for the ratio eB/T^, and we compare the results with the exact ones. 

From the results shown in Tab. HI we see that the EM approximation produces results that already at first order 
have good accuracy compared to the exact values coming from the full Landau sums. It performs particularly very 
well for values of eB ^ T^, reaching convergence within just the very first few terms in the expansion. However, for 
values eB > T^, increasing the order in the EM series, it quickly looses accuracy. This can be traced to the fact that 
as we go to a higher order in the derivatives appearing in the sum on the RHS of Eq. (|4.ip , higher powers of eB /T^ are 
produced, eventually spoiling the approximation (though it oscillates around the true value). This apparent runaway 
behavior is only cured by the introduction of the last term in Eq. (j4.ip . the remainder. 

It should be noted that Eq. (|4.ip is not an approximation to the sum, but it is actually an identity. The important 
term in that context is the remainder, the last term in Eq. (14. ip . In all our numerical tests, within the numerical 
precision for the integrals (for convenience, we have performed all the numerical integrations in Mathematica 43]), 
we have verified that Eq. (|4.ip has agreed with the results from the sum over the Landau levels in the LHS of Eqs. 
(|4.4p and (|4.6p . for all values tested and at all orders (when including the remainder), with a precision of less than 
0.005%. This is quite impressed, recalling that we can in principle keep only the first few terms in Eqs. (|4.4p and 
(|4.6p (the first order approximation) when including the remainder. For numerical computation this is a tremendous 
advantage, since Eq. (|4.ip can easily be coded, compared to having to perform sums over very large Landau levels 
(particularly in the cases of low magnetic fields) and the multiple integrations that may be required. 

Given all the advantages of using the EM formula, it is quite surprising that it is not used frequently in physics, 
despite being well known in mathematics. In particular, we see that it is well suitable in problems involving external 
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TABLE I: The precision of the EM formula (without the remainder) at each order (up to 4th-order), compared with the result 
from the Landau sum. In all cases we have considered Q, — in Eqs. (|4.4|) and (|4.6|) . when extending them up to 4th-order. 



magnetic fields. To our knowledge, one of the few works that we are aware of that have previously made use of 
the EM formula before were for example Refs. [43 - l46j . In Ref. [3l the author derived the effective potential of the 
abelian-Higgs model, including both thermal effects and an external magnetic field, in the one-loop approximation. 
The EM formula was used at its zeroth order (by just transforming the sum in an integral) to obtain analytical 
results for the effective potential in the high magnetic field region. But from the results of Tab. HI this is exactly the 
region and order in the approximation that lead to the largest errors. In (45j the EM formula was used to obtain an 
expression for the effective potential for vector bosons and to study pair production in an external magnetic field, 
while in [4^ it was used to obtain analytical expressions for the internal energies for non-interacting bosons confined 
within a harmonic oscillator potential. None of these previous works have assessed the reliability of the use of the 
EM formula. 

As an application of the EM formula to our problem, we estimate from it the dependence of the critical temperature 
with the magnetic field. Recall from Eq. p. 131) the critical temperature can be obtained by setting 51 = and if = 
in the equation for the VEV. Thus, from Eq. p. 261) . we need to solve the equation 



6— -4X(r„B) 







(4.7) 



As we have seen from the results of Tab. U the EM formula for low magnetic fields produce reliable results already at 
the second order, and for very small magnetic fields, good precision is reached even at first order. Thus, for example 
from the expansion in Eq. (|4.6p . and performing the k integral in the first term, we obtain 
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LX ~ / In ( 1 - e-'^^^+^ 

1 1 1 



2n Jo Vz^ + a2 e^^^ - 1 



eB 1 /I . 

where = eB/T^. Results for the integrals in Eq. (|4.8p can be found in the low magnetic field regime, where a <C 1. 
This regime is equivalent to the high temperature regime in quantum field theory, from where a series expansion in 
a can be found. In fact, the integrals in Eq. (14. 8p ca n be easily related to the hi{a) Bose-Einstein integrals found in 
that context (see for example the App. A in Ref. |14j). from where we obtain the leading order terms in the expansion 
in powers of a: 



Going to higher order in the EM series only lead to 0(a") (with n > 0) terms, for a — \J eB /Tc <^ 1. 

Using the result (|4.9p in Eq. p.23|) and then back in Eq. (|4.7I) . after some algebra, we obtain the approximate 
result for the critical temperature as a function of the magnetic field. 



n, , 25eB ( I 11527r2TO2\ 

- + 32^ + V ^ + ^5A^ j ■ ^'-''^ 

The result (j4.10p shows the growth of the critical temperature with the magnetic field, as seen numerically from the 
results obtained in the previous section. A comparison with our previous results also shows that Eq. ()4.10p provides 
an excellent fit for Tc{B), leading to results with less than 1% error compared to the full numerical results for Tc, for 
values of field such that eB/T^ < 1. For example, for the parameters A = 0.1 and eB = 400m2, the approximation 
(|4.10p gives Tc/m ~ 17.97, while the full numerical calculation gives Tc/m = 18.11. 



V. CONCLUSIONS 



In this paper we have revisited the phase transition problem for the self-interacting complex scalar field model 
when thermal effects and an external magnetic field are present. We have studied this problem in the context of 
the nonperturbative method of the optimized perturbation theory. We have shown that the OPT method preserves 
the Goldstone theorem and that carrying out the approximation to first order is analogous to the ring diagram 
resummation method. The OPT carries out the resummation of self-energy diagrams in a self-consistent way, avoiding 
overcounting issues, which have previously plagued the ring diagram method in its very early applications. 

By using the OPT method, we have demonstrated that the phase transition in the model is always second order in 
the presence of thermal effects and by including an external magnetic field, there is no change to the phase transition 
order. The effect of the external magnetic field is to strengthen the symmetry broken phase, producing a larger VEV 
for the field and also a larger critical temperature. The study we made in this paper considered only the effects of 
the magnetic field and temperature on the charged scalar field (the model is one with a global U{1) symmetry). Our 
study, in a sense, is closer to the case of studying the effects of B and T in the linear sigma model [l^, where no 
couplings to gauge fields (except to the external field) are considered in general. The inclusion of other field degrees 
of freedom, or when promoting the symmetry from a global one to a local symmetry, can of course change both 
qualitatively and quantitatively how the magnetic field affects the system. In particular, as concerning the possibility 
of producing a first order phase transition, either due to thermal effects, or by a magnetic field or from both. This is 
what we expect in the context of the scalar QED, where due to the screening of the magnetic field in the broken phase, 
the phase diagram and transition are very different. The phase structure can also be very different when adding other 
interactions and having other symmetries. For example, in the context of the electroweak phase transition [9l.[l0ll48|. 
it has been shown that the effect of the magnetic field is to make the first order transition stronger. By omitting 
vacuum energy terms from the effective potential, which may have phenomenological motivations, as in Ref. [Zt'I, may 
also lead to a first order transition, instead of a second order one. 
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Still comparing our results with those obtained in the abelian Higgs model, there is a tantalizing question of whether 
for very strong magnetic fields a new phase could be formed. This could be for example a phase with global vortices 
condensation, thus restoring the symmetry again. This would be in analogy to the local, Nielsen- Olesen vortex 
conndensation that can happen in type II superconductors [4(?]. Vortex condensation in that case is energetically 
favorable to happen for fields beyond a critical value and when the mass of the Higgs field becomes larger that the 
mass of the gauge field. This is an interesting possibility to investigate in the future. 

As an aside, but a complementary part of this work, we have verified the reliability of the use of the Euler-Maclaurin 
formula as an approximation for the sum over Landau energies in a magnetic field. We have verified that it produces 
results with errors of less than 0.1% already at the first few orders in the EM formula and that it leads to particularly 
suitable approximations in the low magnetic field region {eB/T"^ <^ 1). The low magnetic field region is typically the 
region where we have to face the problem of summing over very larger number of Landau levels, which, in practice, 
can be overwhelming in terms of CPU time, when working numerically. The EM formula in this case can be a very 
valuable tool, both for numerical computations, but also for obtaining approximate analytical results, which would 
be otherwise very difficult through the use of the Landau sums directly. As an application of the EM formula, we 
have obtained an approximate expression for the dependence of the critical temperature with the magnetic field in 
our problem. We expect that the EM formula can also be used in many other applications involving the effects of 
magnetic fields in phase transitions, like in condensed matter problems or in the recent interest in studying the effects 
of magnetic fields in the QCD phase transition. Furthermore, our results (including the use of the OPT method) can 
also be of interest in the study of cosmological phase transitions in general. Magnetic fields can be easily generated in 
the early universe [1^ . These fields can then influence subsequent cosmological phase transitions or also be important 
in particle physics phenomena in the early universe. Works in those contexts are in progress and we will report on 
them elsewhere. 
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Appendix A: Renormalization in the OPT method 



In this appendix we briefly explain the renormalization of the effective potential Eq. (j2.9l) . First, note that the 
interpolation procedure in the OPT method, Eq. (j2.5l) . introduces only new quadratic terms, thus, it does not alter 
the renormalizability of the original theory. Therefore, the counterterms needed to render the theory finite have the 
same polynomial structure of the original Lagrangian J3Ci] . 

In obtaining the renormalized effective potential, we first note that from the order S term Eq. (j3.4p . we obtain the 
mass renormalization counterterm entering in ALct,s in Eq. (j2.5p . 



Iasi<I>! , (Al) 



where 



From Eq. (|Aip . we obtain the two mass counterterms diagrams contributing to the effective potential at first order 
in the OPT and shown in Fig. [1] 

By collecting all the divergences from the vacuum loop terms in Eq. (|2.9p . we have that the divergence from the 
contribution (|3.4p is canceled by the counterterm Eq. (jAip . A potential temperature and magnetic field dependent 
divergence, proportional to either X in the two- loop vacuum term Eq. (|3.8p . or X in the case of a finite external 
magnetic field, is explicitly canceled against identical term coming from the mass counterterm diagrams, Eq. (13.61) . 
The remaining divergences in the effective potential are only vacuum ones, independent of the background field, and 
they can be all canceled by introducing a vacuum counterterm A Vet, added to the effective potential. At the first 
order in the OPT, AV^t is given by 
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(A3) 



Thus, at first order in the OPT, we only require the two counterterms, Eqs. (|Aip and (|A3[) . Going to second order 
in the OPT it will also require a coupling constant renormalization counterterm, Cs2Xipf from which we can built 
vertex counterterm vacuum diagram contributions to the effective potential [33| (see also the first two references in 
(30j). At second order in the OPT we would also obtain 6^ contributions for Eqs. (lAip and (IA3p . Going to higher 
orders in the OPT only produces additional 0((5"), n > 2, contributions to the mass, vertex and vacuum counterterms. 
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